Spin wave dispersion in La2Cu04 
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We calculate the antiferromagnetic spin-wave dispersion in a half-filled Hubbard model for a 
two-dimensional square lattice, and find it to be in excellent agreement with recent high-resolution 
inelastic neutron scattering performed on La2Cu04 [Phys. Rev. Lett. 86, 5377 (2001)]. 
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After over a decade of intense research on the micro- 
scopic origin of high-temperature superconductivity in 
cuprates, there is no general consensus on the microscopic 
Hamiltonian suitable for describing these materials. Nev- 
ertheless, it appears that magnetic fluctuations must play 
an important role. Therefore, the study of magnetic fluc- 
tuations in the high-temperature superconductor parent 
compounds, such as La2Cu04, is an important field of 
research, both theoretical and experimental. 

In two recent papersj^l'H high-resolution inelastic 
neutron-scattering measurements were performed on two 
different two-dimensional spin-1/2 quantum antifcrro- 
magnets. These are copper deuteroformate tetradeuter- 
ate (CFTD) and La2Cu04. Surprisingly, the dispersion 
at the zone boundary observed in the two materials, 
does not agree with spin- wave theory predictions J3 More- 
over the amount of dispersion is not the same for both 
materials. In CFTD the dispersion is about 6% from 
a;(7r/2,7r/2) to aj(7r, 0), whereas in La2Cu04 it is about 
-13% along the same direction. In the case of CFTD the 
dispersion at the zone boundary can be explained using 
the nearest-neighbor Heisenberg model alone,a and high- 
precision quantum Monte Carlo simulations have con- 
firmed that this is sod On the other hand, an explanatioH. 
for the observed dispersion in La2CU|Q4 was proposed,EJ 
using an extended Heisenberg modelQu involving first-, 
second-, and third-nearest-neighbor interactions as well 
as interactions among four spins. This extended model 
was obtained from the Hubbard model, using perturba- 
tion theory, and is diagonalized afterward using classical 
(large-S*) linear spin-wave theoryJH 

The La2Cu04 results clearly show that the usual 
Heisenberg model is insufficient to explain the experi- 
mental data, and that the Hubbard model is the correct 
Hamiltonian |£or describing the magnetic interactions in 
the cupratesfl In this work we do not use perturbation 
theory for deriving an effective magnetic Hamiltonian. 
Instead we work directly with the Hubbard model. We 
consider a half-filled Hubbard model in a spin-density- 
wave (SDW)- broken symmetry ground state and, by 
summing up all ladder diagrams, we compute the trans- 
verse spin susceptibility and from this obtain the spin- 
wave dispersion. 

The Hubbard model for a square lattice of N sites is 



defined as 
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where the sum over counts each pair of nearest 

neighbors only once, the momentum k runs over the 
Brillouin zone, the electronic energy dispersion e(fc) — 
~2tcoskx ~ 2t cos ky has the nesting vector Q = (7r,7r), 
and 

k,k' ,q 

The broken-symmetry state is introduced by consid- 
ering the existence of an off-diagonal Green's function 
given by 

FM r-r')^ -{TrC^^±Qjr)cljT')} . (3) 
in addition to the usual Green's function: 

G.(p; T - r') - -(r.cpv(r)4,,(r')) . (4) 

In mean-field theory for a SDW, the two propagators 
are related as shown in the Feynnian diagrams depicted 
in Figs. 1^ and |^. Single (doubled) arrowed lines repre- 
sent diagonal (off-diagonal) propagators. The single line 
represents the free propagator, whereas double lines rep- 
resent mean-field propagators. The Hubbard interaction 
is represented by a dashed line. 
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FIG. 1. Feynman diagrams corresponding to the SDW 
mean-field equation for Ga{p', t — t). 
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FIG. 2. Feynman diagrams corresponding to the SDW 
mean-field equation for Fcr(p; t ~ t') 

The solution to the mean-field equations yields the 
SDW staggered magnetic moment m, which is defined 
as 
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The staggered magnetic moment m is reduced from its 
Neel value for finite values of t/U, and its behavior as 
function of t/U at zero temperature is shown in Fig. ^ 




FIG. 3. Staggered magnetization as function of U /t at 

r = o K. 

Both the above Green's functions have two-pole struc- 
tures given by 



iujn - {p) iun - {p} 



(6) 
(7) 
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where E{p) = - ^(P + QW + U^m^, and^(p) = 

e{k) — ji^ and the coherence factors read 



iun - E+ {p) iujn ~ E- (p) 
where the energies E± are given by 

E±{p) = + ^ 2 ''^ ' 
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It is known that the mean-field treatment of the spin 
dynamics of itinerant strongly correlated electronic sys- 
tems yields satisfactory results, as in the study of the 
dynamic spin response function of the cuprates' super- 
conducting state.u In order to describe the spin dynamics 
of the system, we consider the transverse spin suscepti- 
bility X 'i^n), which is defined as 
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where f3 = l/T is the inverse temperature, Tr is 
the chronological order operator (in imaginary time), 

S-{q) = Ep4iCj;+?,h and S+iq) = [S' {q)]K The 
above expression can be written as 

n^a^o Jo 
where d.c. stands for differently connected diagrams. 
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FIG. 4. The four first-order bubble diagrams for 
x\ (<f, iWn) in the SDW state. A summation over the in- 
ternal momentum along the interaction line is implied. 
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We now compute x *'^n) by summing up all lad- 
der diagrams and taking into account the existence of 
two possible Green's-function lines. The first-order bub- 
ble diagrams are shown in Fig. |4[ The complete ladder 
summation is given by 



400 



ladder 



{q, iu!„ 



-jx + y)[l + Xjx + y)] + A(zi + za)^ 
[1 + X{x + y)][l + Xix + y)] - X^z^ + Z2Y ' 

(13) 



with A = U/N and 



x{q, iw„) = i ^ Gj, (p, iun ')G| [p + q, ii^n ' + jw„) , 
y{q, iujn) = -p^Fi [p, iuJn ')Fi [p + g, iw„ ' -|- iw„) , 

zi {q, iw„) = i ^ Gx (p, iU!n ')^T (P + 9' 'i'^n ' + iuJn) , 
n ' ,p 

Z2(g, iWn) = ^ X! ')'^T (P + 9, ' + i^n) , 

x{q, itJ„) = i ^ Gx (p, iw„ (p + g + Q, ' + it^„) . 



The retarded susceptibility is obtained from Eq. ( |13D , 
performing the analytical continuation iw„ ^ + iO+. 
The poles of the retarded susceptibility give the energy 
ioiq) of the spin excitations of the system as well as their 
lifetimes. The equation for the poles is 



{\^X\xiq,uj)+y{q,Uj)W{\^X\x{q,uj)^y%uj)\} 
-X\zM^)^zMuj)f (14) 



It is not possible to solve Eq. ( [Til ) analytically for arbi- 
trary \J and t. In the limit tjU ^ at half-filling, from 
( p^ ) we recover the same result as in linear spin-wave 
theory for the nearest-neighbor Heisenberg model, 



j 

i^iqx,qy) ;7^V^ ~ (cosg^r + COS qy)"^ , (15) 



which predicts that a;(0,7r) = aj(3r/2, 7r/2), in disagree- 
ment with the experimental datap and m ^ 1, showing 
no magnetic moment reduction from the Neel state value. 
Of course Eq. (^ holds asymptotically for t/U <^ 1. We 
remark that the factor 1 / 77i,plays the role of the quantum 
renormalization factor Zr.u 




FIG. 5. Spin-wave dispersion, in meV, 
along high-symmetry directions in jthe Brillouin zone. The 
circles are the data reported in RefJj at 10 K. The solid line 
is the analytical result (at K) for U = 1.8 eV and t — 0.295 
eV. The momentum is in units of 2n, and M — (1/2, 1/2), 
X = (1/2,0), and T = (0,0). 



On the other hand, for finite t/U, Eq. (^^ has to 
be solved numerically. Considering the half-filled case 
(p = 0), appropriate for La2Cu04, we computed the 
spin-wave dispersion u!{q) for U = 1.8 eV and t = 0.295 
eV along the high symmetry directions, in the two- 
dimensional Brillouin zone, considered in Ref. 0. These 
values agree with those used in Ref. U = 2.2 ± 0.4 eV 
and t — 0.30±0.02 eV. The results (solid line) are given in 
Fig. H together with the experimental results (in circles) 
at r=:10 K. It is clear that t^(0,7r) > uj{TT/2,n/2), that 
is, a dispersion at the zone boundary is obtained. For 
these values of t and U the staggered magnetic moment 
is m = 0.832, and therefore 1/m — 1.20, which agrees 
well with Zc = 1.18 used to fit the data in Ref. |l|. 

In their interpretation of the experimental data, the 
authors of Ref. |l] considered an effective Hamiltonian in- 
corporating ring exchange. In such a model the electron 
not only makes a virtual trip to its nearest neighbor, but 
also makes virtual excursions around a loop visiting its 
second neighbors. If we had written the transverse sus- 
ceptibility in coordinate space, it would be clear that in 
a ladder summation the electron goes around larger and 
larger rings before it comes back to the original site with 
its spin dipped. Therefore, such a good agreement be- 
tween the perturbation theory calculation in Ref. |^ and 
our ladder summation is not surprising. Althought we 
have not presented the results here, Eq. ( 14) also predicts 
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a continuum of excitations which are gapped relatively to 
the ground state. 

In summary, we have shown that a ladder summation 
based on the SDW state can account satisfactorily for 
the measured spin- wave dispersion of La2Cu04 at all en- 
ergies. The quality of the fitting points out that it is not 
needed to derive an effective spin Hamiltonian from the 
Hubbard model in order to obtain agreement with the 
data. 

We gratefully acknowledge G. Aeppli for discussions, 
and for providing us the experimental data we show in 
Fig. ||. We thank A. Sandvik for calling our attention to 
this problem. 



1 R. Coldea, S. M. Hayden, G. Aeppli, T. G. Perring, C. D 
Frost, T. E. Mason, S.-W. Cheong, and Z. Fisk, Phys. Rev 
Lett. 86, 5377 (2001). 

^ H. M. R0nnow, D. F. McMorrow, R. Coldea, A. Harrison 
I. D. Youngson, T. G. Perring, G. Aeppli, O. Syljuasen, K 
Lefmann, and C. Rischel, Phys. Rev. Lett. 87, 37202 (2001) 

^ J. Igarashi, Phys. Rev. B 46, 10763 (1992). 

^ Anders Sandvik and Rajiv R. P. Sing, Phys. Rev. Lett. 86 
528 (2001). 

^5 M. Takahashi, J. Phys. C 10, 1289 (1977). 

A. H. MacDonald, S. M. Girvin, and D. Yoshioka, Phys 

Rev. B 41, 2565 (1990); 37, 9753 (1988). 
^P. W. Anderson, Science 255, 1196 (1987). 



' See, for example, Jan Brinckmann and Patrick A. Lee, cond 



mat/0110316 [J. Low Temp. Physics (to be published)], and 



references therein. 



4 



